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I.  INTRODUCTION 


The  dependence  of  the  threshold  voltage,  subthreshold  I-V 
characteristics,  and  radiation  response  of  n-channel  AlGaAs/GaAs  MODEFTs 
(modulation  doped  field-effect  transistors)  on  acceptor  doping  density  has 
been  described  previously  (Refs.  1 -A ) .  These  analyses  have  now  been 
extended  to  describe  the  dependence  of  MODFET  high  channel  carrier  density 
I-V  characteristics  on  acceptor'  doping  density.  The  effect  of  acceptor 
doping  density  on  the  experimental  threshold  voltage  and  device  capacitance 
is  explicitly  addressed. 

The  band  structure  of  a  typical  AlGaAs(n)/GaAs  hetero junction  with 
Schottky  barrier,  $  ,  at  the  gate,  and  a  spacer  layer  at  the  interface 
under  bias  V  ,  is  shown  in  Fig.  1.  In  the  depletion  layer  approximation, 

D 

the  donors  and  acceptors  are  assumed  to  be  completely  ionized  in  the  doped 

AlGaAs  layer  d,  the  spacer  layer  a,  and  in  the  depletion  layer  W.  The 

doping  densities  ND  and  are  assumed  constant.  The  quasi-two-dimensional 

electron  eigenstates  at  the  interface  are  solved  for  using  a  trian'  jlar 

potential  well,  and  only  the  lowest  subband  is  included  in  the  calculation. 

A  delta-function  channel  charge  distribution  at  the  average  channel  width 

is  assumed.  Band  bending  from  the  interface  at  (d  +  a)  to  the  edge  of  the 

depletion  region  (W  +  d  +  a)  is  the  difference  of  the  position  of  the 

conduction  band  relative  to  the  Fermi  level  in  the  GaAs  far  from  the 

interface  (E_/2  +  <t>,  , .  ),  and  the  Fermi  level  Er  relative  to  the  bottom  of 

g  bulk  ’  i 

the  two-dimensional  channel. 

Under  the  restrictions  imposed  by  these  assumptions,  Poisson's 
equation  may  be  integrated  across  the  structure  to  obtain  the  applied  gate 
voltage  as  a  function  of  device  geometry,  doping  densities,  and  channel 
charge  ns: 

V  =  *  -  AE  ♦  Er  +  (q/e)[-N  d2/2  +  (N  W  +  n  )(d  +  a)]  (1) 

m  o  i  u  h  s 
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where 


AE  -  conduction  band  offset 

In  the  one-subbarid  approximation,  the  Fermi  level  Ep  may  be  calculated  as 

Ep  -  Ejj/q  +  (  kT.'q )ln[  exp(  uh^n3 - m^kT )  -  1]  (2) 

where  the  first  quantum  level  Eq  may  be  calculated  in  the  tr iangui ar-wel 1 
approximation  and  is  given  by 

E0  =  (  9+<2/8m1 )  (MEjntp  'tf2)2/3  (3) 

The  field  E(  at  the  interface  is 

Ex  =  (q/e)(ns  -  Nfl )  (4) 

Substitution  of  Eqs.  (2-4)  into  Eq .  (1)  yields 

vg  =  v0  *  f(ns)  (5) 

where 

Vn  =  $  -  AE  - ( q  ' *  )Nnd2/2  (6) 

0  m  c  L) 

The  function  f(ns)  may  be  written  as 

f(ns)  -  ( q/e) (d  ♦  a)(NftW  «•  ns)  +  CQ(NAW  ♦  ns)2'3 

+  (kT/q)ln[exp(ns/nc)  -  1]  (7) 
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where 


Cq  =  (9tf^/8rrijq)  (4q2mj/tt~e  P/3 

=  [-1.7  x  10~9  V-cn/473] 

/K"  -  Planck  constant  divided  by  2n 

m ^  =  longitudinal  effective  mass  of  the  carriers 

We  have  retained  terms  containing  the  depletion  width  W  in  the  definition 
of  f(ns)  because  the  depletion  width  is  implicitly  a  function  of  the 
channel  charge  n  . 

The  quantities  q,  e,  k,  and  T  are  the  elemental  charge,  AlGaAs(GaAs) 
permittivity  (assumed  identical),  Boltzmann  constant,  and  absolute 

temperature.  * 

Similarily,  charge  density  nc  is  a  function  of  physical  constants  and 
the  effective  carrier  mass:  » 


nc  =  n^/mjkT  (9) 

and  is  equal  to  -8. A  *  10  cm  . 

A  discussion  of  the  dependence  of  the  depletion  widtn  W  on  acceptor 
density  has  been  given  elsewhere  (Ref.  1).  In  section  II  we  describe  the 
mathematical  properties  of  the  function  f(ns),  in  the  high  channel  carrier 
density  region,  and  exploit  the  results  to  describe  the  characteristics  of 
these  devices. 


» 
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II.  HIGH  DENSITY  REGION 


A.  DEFINITION 

We  define  the  high  density  region  such  that  ng  >  nc  over  the  whole 

channel.  In  this  region  Eq.  (7)  may  be  expanded  in  a  Taylor  series  in  ns 

about  nc.  The  results  of  this  approximation  to  first  order  are  shown  in 

Fig.  2  for  two  extremes  of  acceptor  doping  density.  The  solid  lines  are 

the  results  of  Eq.  (7),  and  the  dashed  lines  are  the  result  of  the  first 

11  P 

order  expansion  in  ns  about  n  .  Above  nc  (8.4  x  10  cm  ) ,  the  expansion 

11  2 

is  quite  good.  Much  below  nc  (<  3.0  x  10  cm  ) ,  the  first  order 

expansion  departs  from  the  exact  result  and  approaches  a  constant. 

Substitution  of  the  first  order  expansion  for  f(r.s)  in  Eq.  (5)  and 

inverting  to  find  n  as  a  function  of  V  yields 

^  6 

ns  --  nc  ♦  -  vo  -  f(nc)l/(kT/q)  (10) 

i 

where  K  is  a  constant  that  depends  on  the  device  geometry,  doping 
densities,  depletion  width,  and  physical  constants.  This  form  for  ns  is 
different  than  previously  assumed  (Ref.  5),  which  ignores  the  contribution 
from  n£  and  from  f(nc)  and  assumes  that  the  reciprocal  of  K  is  the  AlGaAs 
layer  capacitance  per  unit  area  per  unit  charge  times  the  absolute 
temperature  in  electron  volts.  In  our  formulation,  near  saturation,  K  is 
given  by: 

K  =  ( q/kT ) ( q/e ) ( d  +  a)  ♦  (2/3)C0(q/kT)/(NftW  +  np)1/3  +  1.58/nc  (11) 

B.  DEVICE  CAPACITANCE 

The  derivative  of  Eq.  (10)  with  respect  tc  V  yields  the  device 

•l  O 

capacitance  per  unit  area,  which  may  be  written  in  the  following  form 
(Ref.  6): 

i 

Carea  =  c/(d  +  a  +  Ad>  (  12) 
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Function  f(ns)  versus  Channel  Carrier  Density  for 
Two  Extremes  of  Acceptor  Density.  The  solid  curves 
are  the  exact  results  of  Eq.  (2);  the  dashed  lines 
are  the  results  of  the  linear  approximation,  described 
in  text. 


where  Ad  i a  invoked  to  account  for  the  discrepancy  between  the  AlGaAs  layer 
capacitance  e'(d  +  a)  and  the  device  capacitance.  Taking  the  derivative  of 
Eq .  (10)  and  using  Eq.  (12)  to  solve  for  Ad  yields 

Ad  =  (  Al  .  3q)Cn  ( N„W  +  n_  ) 1  3  +  1 .58(  kT 'q ) ( t ' q )  n  (13) 

For  low  acceptor  densities  (■  101-*  err,"3),  Ad  levels  off  at  a  value  of  -89  A, 
which  is  in  good  agreement  with  values  cited  in  the  literature  (Ref.  5). 

As  the  acceptor  density  increases,  Ad  decreases.  At  an  acceptor  density  of 
1017  cm"-5,  Ad  decreases  to  -74  A. 

C.  EXPERIMENTAL  THRESHOLD  VOLTAGE 

MODFET  threshold  voltages  are  determined  experimentally  by  extrapolat¬ 
ing  the  saturation  current,  or  square  root  of  the  saturation  current, 
versus  gate  voltage  to  zero.  The  gate  voltage  intercept  is  the  experiment¬ 
ally  determined  threshold  voltage.  We  may  approximate  the  experimental 
threshold  voltage  by  solving  Eq.  (10)  for  the  gate  voltage  when  n  is  equal 
to  zero.  This  yields  a  threshold  voltage  which  differs  from  the  strong 
inversion  definition  of  threshold  voltage  (Ref.  1).  The  difference  between 
this  approximation  for  the  threshold  voltage  and  the  strong  inversion 
threshold  voltage  is  given  by 

»th  =  r<"c>  -  -  «"c  carea  <'«> 

where  the  channel  charge  at  threshold  nt(l  is  equal  to  the  acceptor  density 
N^ ,  times  the  average  channel  width  zav,  which  may  be  calculated  in  the 
triangular-well  approximation  using  variational  functions  (Refs.  7,  8). 

This  definition  of  threshold  has  been  described  elsewhere  (Ref.  1). 

In  Fig.  3  the  threshold  voltage  difference,  Eq.  (14),  is  plotted 

versus  acceptor  density.  Each  term  in  Eq.  (14)  is  plotted  separately.  At 

14-7 

low  acceptor  densities  (  ■•  10  cm  -*)  the  difterence  may  be  as  much  as 
0.25  V.  This  difference  decreases  as  the  acceptor  density  increases.  As 
this  difference  depends  on  the  acceptor  density,  a  comparison  of 


experimental  threshold  voltages  may  not  be  appropriate  if  acceptor 
densities  differ  by  a  significant  amount. 

D.  I-V  CHARACTERISTICS 

In  the  gradual  channel  approximation,  charge  control  is  determined  by 
the  effective  potential  in  the  channel: 

V(x)  r  Vg  -  V c ( x )  (  15) 

where  Vc(x)  is  the  channel  voltage  under  the  gate  at  point  x.  Using  Eq. 
(10)  we  solve  for  the  carrier  density  in  the  channel: 

ns(x)  =  nc  +  K_1[Vg  -  VQ  -  f(nc)  -  Vc(x)]/( kT/q)  (16) 

The  form  of  Eq .  (16)  allows  the  source-drain  current  to  be  calculated  in 
the  usual  way  (Ref.  5).  The  result  is 

I HD  =  lU/L)u[lnc  «  l(Carea/q)(Vg-V0-r(nc)]||Vc(L,-Vc(0)l 

-  <Carea/‘l)lUc(L)-'lc<0>)l  (17> 

where  u  is  the  channel  mobility  and  Z/L  is  the  gate-width  to  gate-lengch 
ratio.  For  a  grounded  source,  in  the  limit  of  zero  source-drain  resis¬ 
tance,  we  recover  the  usual  dependence  on  drain  voltage. 

In  the  development  of  Eq.  (17)  we  have  assumed  that  each  point  in  the 
channel  has  a  carrier  density  greater  than  n  .  This  places  a  limit  on  the 
bias  conditions  for  the  applicability  of  E}.  (17).  Evaluating  Eq.  (16)  at 
the  drain  contact  and  using  the  condition  that  n  (L)  must  be  greater  than 
nc ,  yields  the  following  limitations  on  the  bias  conditions  for  a  grounded 
source  and  zero  source-drain  resistance: 


(YV 


f  (n. 


(  18) 
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where  Vp,  is  the  drain  voltage.  Shown  in  Fig.  4(a)  is  I^g  versus  Vp,  for  a 
grounded  source  and  zero  source-drain  resistance  for  various  acceptor 
doping  densities  and  the  device  parameters  shown.  The  heavy  dots  indicate 
the  point  at  which  the  calculation  is  no  longer  valid  by  virtue  of  Eq.  (18). 
As  the  acceptor'  density  increases,  the  current  at  a  given  drain  voltage 
decreases  and  the  limit  of  the  appi i cabi l i ty  of  the  calculation  decreases, 
in  drain  voltage,  for-  a  given  gate  voltage.  A  similar-  graph  is  shown  m 
Fig.  4(b)  in  which  lj.^  versus  Vp,  is  shown  for  a  given  acceptor  density  and 
various  gate  voltages.  As  the  gate  voltage  decreases,  the  current  at  a 
given  drain  voltage  decreases,  as  expected,  and  th>  egiori  oi  applicability 
of  the  calculation,  in  drain  voltage,  decreases,  inis  result  implies  that 
using  an  equation  of  the  form  given  in  Eq.  (17),  which  results  from  a 
linearization  of  f(n  )  above  n  ,  to  infer  the  nature  of  the  saturation 
characteristics  in  MODFETs ,  in  which  the  drain  portion  of  the  channel  has  a 
channel  carrier  density  much  less  than  n  ,  is  suspect  (Ref.  5). 
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DRAIN  VOLTAGE  [Volts] 


H  -  4500  cm‘/V  s 
Nd  =  1  x  10’8  cm  3 
NA  -  1  x  10'"  cm  ' 


DRAIN  VOLTAGE  [Volts) 


Drain-Source  Current  versus  Drain  Voltage  for 
Various  Acceptor  Densities  (a),  and  Various 
Gate  Voltages  (b).  Heavy  dots  indicate  the 
limits  of  the  calculation  as  given  by  Eq.  (11) 
in  the  text. 


HI.  SUMMARY 


I 

We  have  developed  a  triangular-well,  one-subband  depletion  layer  modci 
to  describe  the  high  channel  density  operation  of  MODFETs.  Trie  effects  of 
acceptor  density  on  the  I-V  characteristics  in  the  high  channel  density 
region  has  been  investigated.  The  depletion  layer  charge,  due  to  ionized 
acceptors,  is  shown  to  account  for'  the  discrepancy  between  the  device 
capacitance  and  the  Albans  layer-  capac  i  tehee ,  as  described  in  trie  litera¬ 
ture  (Ref.  6).  The  depletion  layer  charge  also  accounts  for  the  discrep¬ 
ancy  between  the  theoretical  thresnoid  voitage,  in  the  strong  inversion 
model  (Ref.  1),  and  the  experimental  threshold  voltage,  as  determined  by 
extrapolation  of  the  source-drain  current.  Therefore,  comparison  of 
threshold  voltage  characteristics,  particularly  for  devices  with  appreci¬ 
ably  different  acceptor-  densities,  should  be  based  on  a  consistent 

* 

description  of  threshold  that  accounts  for  ionized  acceptors,  i.e.,  the 
strong  inversion  model. 

t 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer”  for  national  security 
projects,  specializing  in  advanced  military  space  systems.  Providing  research  support,  the 
corporation’s  Laboratory  Operations  conducts  experimental  and  theoretical  investigations  that 
focus  on  the  application  of  scientific  and  technical  advances  to  such  systems.  Vital  to  the  success 
of  these  investigations  is  the  technical  staffs  wide-ranging  expertise  and  its  ability  to  stay  current 
with  new  developments.  This  expertise  is  enhanced  by  a  research  program  aimed  at  dealing  with 
the  many  problems  associated  with  rapidly  evolving  space  systems.  Contributing  their  capabilities 
to  the  research  effort  are  these  individual  laboratories: 

Aerophysics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat  transfer 
and  flight  dynamics;  chemical  and  electric  propulsion,  propellant  chemistry,  chemical 
dynamics,  environmental  chemistry,  trace  detection,  spacecraft  structural  mechanics, 
contamination,  thermal  and  structural  control;  high  temperature  thermomechanics,  gas 
kinetics  and  radiation;  cw  and  pulsed  chemical  and  excimer  laser  development, 
including  chemical  kinetics,  spectroscopy,  optical  resonators,  beam  control,  atmos¬ 
pheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmospheric 
optics,  light  scattering,  state-specific  chemical  reactions  and  radiative  signatures  of 
missile  plumes,  sensor  out-of-field-of-view  rejection,  applied  laser  spectroscopy,  laser 
chemistry,  laser  optoelectronics,  solar  cell  physics,  battery  electrochemistry,  space 
vacuum  and  radiation  effects  on  materials,  lubrication  and  surface  phenomena, 
thermionic  emission,  photosensitive  materials  and  detectors,  atomic  frequency  stand¬ 
ards,  and  environmental  chemistry. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device  physics, 
compound  semiconductors,  radiation  hardening;  electro-optics,  quantum  electronics, 
solid-state  lasers,  optical  propagation  and  communications;  microwave  semiconductor 
devices,  microwave/millimeter  wave  measurements,  diagnostics  and  radiometry,  micro¬ 
wave/millimeter  wave  thermionic  devices;  atomic  time  and  frequency  standards; 
antennas,  rf  systems,  electromagnetic  propagation  phenomena,  space  communication 
systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals,  alloys, 
ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  nondestructive 
evaluation,  component  failure  analysis  and  reliability;  fracture  mechanics  and  stress 
corrosion;  analysis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures 
as  well  as  in  space  and  enemy-induced  environments. 

Space  Sciences  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray  physics, 
wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and  ionospheric 
physics,  density  and  composition  of  the  upper  atmosphere,  remote  sensing  using 
atmospheric  radiation;  solar  physics,  infrared  astronomy,  infrared  signature  analysis; 
effects  of  solar  activity,  magnetic  storms  and  nuclear  explosions  on  the  earth’s 
atmosphere,  ionosphere  and  magnetosphere;  effects  of  electromagnetic  and  particulate 
radiations  on  space  systems;  space  instrumentation. 


